Dressed Feshbach molecules in the BEC-BCS crossover 
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We present the RPA theory of the BEC-BCS crossover in an atomic Fermi gas near a Feshbach 
resonance that includes the relevant two-body atomic physics exactly. This allows us to determine 
the probability Z for the dressed molecules in the Bose-Einstein condensate to be in the closed 
channel of the Feshbach resonance and to compare with the recent experiments of Partridge et al. 
||cond-mat/0505353 with ^Li. We determine for this extremely broad resonance also the condensate 
density of the dressed molecules throughout the BEC-BCS crossover. 



Introduction. — The superfluid phase in an atomic 
Fermi gas near a Feshbach resonance realizes a funda- 
mentally new state of matter, which shows a macroscopic 
coherence between atom pairs and molecules that is con- 
trolled by the applied magnetic field. As a result, such 
a gas offers the exciting possibility to study in great de- 
tail the crossover between the Bose-Einstein condensa- 
tion (BEC) of diatomic molecules and the Bose-Einstein 
condensation of atomic Cooper pairs, i.e. .the Bardeen- 
Cooper-Schrieffer (BCS) transition In fact, 

the BEC-BCS crossover region is presently already ac- 
tively being explored by a number of experimental groups 
around the world 0, lE 01111 E3- 

In more detail the physics of the BEC-BCS crossover 
occurring near a Feshbach resonance can be understood 
as follows: The superfluid phase of the gas is always asso- 
ciated with a Bose-Einstein condensate of pairs, but the 
wave function of the pairs or dressed molecules is given 
by the linear superposition 

(rlxdrosscd) = VZxm(r) I closed) -|- Vl - Zxaa(r)|open) . 

In the BEC limit the applied magnetic field is taken such 
that the bare molecular energy level lies far below the 
threshold of the two- atom continuum and we have Z ~ 1. 
In that case we are dealing with a Bose-Einstein conden- 
sate of tightly-bound diatomic molecules and the spatial 
part of the pair wave function is equal to the bare molec- 
ular wave function Xm(r)- The spin part of the pair wave 
function is then equal to [closed), i.e., the spin state of 
the closed channel that causes the Feshbach resonance 
In the BCS limit the bare molecular energy level 
lies far above the threshold of the two-atom continuum 
and can be adiabatically eliminated. We then have that 
Z ~ and the spatial part of the pair wave function 
equals the usual BCS wave function for atomic Cooper 
pairs Xaa(r). This Cooper-pair wave function depends 
also on the magnetic field, because the effective attrac- 
tion between the atoms after the adiabatic elimination of 
the bare molecular state depends on the energy of that 
state. The spin state of the Cooper pairs is, however, 
always equal to the spin state jopen) of the open channel 
of the Feshbach problem. 

The probability Z plays therefore a crucial role in the 



description of the BEC-BCS crossover since it quantifies 
the amount of coherence between the atom pairs and the 
bare molecules in thegas. Unfortunately, however, the 
various theories Ill[ll[ll[ll[ll[ll[ll[llll3 that 
are presently being used to understand the outcome of 
the experiments are unable to accurately determine this 
quantity. This comes about because, either Z is assumed 
to be zero from the outset, the many-body theory does 
not incorporate the two-body Feshbach physics exactly, 
or the theory is able to determine only the total num- 
ber of bare molecules in the gas and thus requires an 
additional assumption about the total number of dressed 
molecules to extract Z. This situation is particularly 
unsettling because of the recent ^Li experiments of Par- 
tridge et al. 0, which have used the photodissociation 
rate to measure the value of Z throughout the crossover 
regime. In view of the above situation it is pressing to 
develop an ah initio many-body theory for the calcula- 
tion of Z . How that may be achieved is the main topic 
of this Letter. 

BEC-BCS crossover theory. — Introducing creation 
and annihilation operators for the bare molecules and 
atoms, the effective grand-canonical hamiltonian of the 
gas with chemical potential becomes 1211 12^ 

77 = y dx V^t^(x) +5-2i^ ^,n(x) (1) 

+ g j dx (x)V'T (x)'!/'! (x) -I- ^|(x)V'|(x)?/:in(x) j 

+ j dx V;|(x)^|(x)^TWV'iW : 

where the two hyperfine states of the atoms are de- 
noted by I ]) and | J,), and the magnetic-moment dif- 
ference A/iinag between the hyperfine states | closed) and 
jopen) = (I ti) — I i]))/^ gives the so-called detun- 
ing from resonance 5 — A/i,nag(73 — Bq). Note that 
the atom-molecule coupling constant g and the back- 
ground scattering length Obg depend on the magnetic 
field B in such a manner that the total scattering length 
a = flbg — mg^ /Anh'^S agrees with the Feshbach reso- 
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nance of interest. This is especially important for the 
broad Feshbach resonance near 834 Gauss that is used in 
all ^Li experiments up to date 

From now on we consider only the most interesting 
region close to resonance, where Obg ^ a and the effec- 
tive interaction between the atoms is dominated by the 
resonant part —g'^/S. This suggests that the last term 
in the right-hand side of Eq. can be neglected alto- 
gether. This is, however, not true because we can neglect 
this term only after we have included its effect on the 
atom-molecule coupling |2^. Physically, the reason for 
this subtlety is that the above Hamiltonian is an effective 
Hamiltonian that is only valid for low energies. However, 
to properly account for the two-body physics near a Fes- 
hbach resonance also high-energy states are needed. The 
main effect of these high-energy states is to renormal- 
ize the atom-molecule coupling to g{)<.) = g/{l + ifcobg), 
where k is the relative momentum of the atoms involved 
in the coupling [lllE5l |. Only after having performed this 
substitution are we allowed to neglect the background in- 
teraction between the atoms. 

I 



Without the background interaction the atomic part of 
the Hamiltonian is quadratic. Using standard functional 
methods the fermionic fields can thus be integrated out 
exactly. This leads to an effective action for the molecules 
that at sufficiently low temperatures has a minimum at 
a nonzero value of {ipm{x)) = y/Zn^, where we intro- 
duced the dressed molecular condensate density nmc- Ne- 
glecting fluctuations at this point leads to a mean-field 
theory of the BEC-BCS crossover. As mentioned in the 
introduction, however, this mean-field theory is unable 
to calculate the probability Z since it only determines 
the bare molecular condensate density . We 

therefore also consider the quadratic fiuctuations around 
the minimum of the effective molecular action, i.e., we 
consider the Bogoliubov theory of the bare molecules. 

As expected with a spontaneously broken U{1) sym- 
metry associated with the presence of a Bose-Einstein 
condensate, the gaussian fiuctuations are determined by 
normal and anomalous self energies of the bare molecules, 
which at zero temperature reduce to 



J {2^f^^^ \inw - hLOa{^'+) - huJa{'k'_) zfiw + fiWa (k' + ) + fiWa (k'_ ) 2e(k') J ' 



Here we have also introduced the BCS dispersion 
tujja.{k) — yj (e(k) — + |g(k)pZrti„c, the bare atomic 
dispersion e(k) — tfk'^/2m, the usual BCS coherence fac- 
tors Ma(k) and Wa(k), and the notation k'j_ = k/2 ± k'. 
In terms of the above self energies the minimum of the 
effective action is determined by the exact Hugenholtz- 
Pines relation 2iJ. = S + ?iSii(0, 0) - ;iEi2(0,0), which 
turns out to be equal to a modified BCS gap equation 
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2hLu{k) 2e(k) 



(3) 



Finally, we also need the equation of state, which we 
for consistency reasons llj obtain by differentiating the 
thermodynamic potential with respect to the chemical 
potential. Including the effect of the fluctuations we ob- 
tain for the total density of atoms 



n = -Tr[Ga] + 2Zn,„c - Tr[G,n] + ^Tr 



G, 



dhY. 



(4) 



where Ga and Gm are the Nambu (2 x 2)-matrix Green's 
functions of the bare atoms and molecules, respectively. 
For a given density and magnetic field the last two equa- 



tions determine only the bare molecular condensate den- 
sity and the chemical potential. Hence, we need a third 
equation to determine also Z. 

Dressed molecules. — Before we derive this missing 
equation, let us first discuss in some more detail the 
physics behind the maybe somewhat unexpected equa- 
tion of state in Eq. Q). The first two terms represent 
the mean-field theory without fiuctuations that is most 
often used in the recent literature 

UlElIIlIIi Be- 
cause of the absence of fiuctuations all the molecules are 
Bose-Einstein condensed and there is no depletion. The 
third term precisely corresponds to this depletion. Fi- 
nally, the fourth term physically describes the dressing 
of the bare atoms and molecules. This can be made 
more clear by reformulation the equation of state in terms 
of dressed atoms and molecules, instead of bare atoms 
and molecules. Since every dressed molecule contains 
two atoms, we expect the contribution 2ni„c from the 
condensate of dressed molecules. Indeed, in the BEC 
limit it can be shown explicitly that the atomic den- 
sity -Tr[Ga] = 2 /dk|t;a(k)|V(27r)3 contains exactly the 
expected contribution 2(1 — Z)ninc of paired atoms in 
the Bose-Einstein condensate of dressed molecules. The 
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FIG. 1: Spectral functions of the bare ®Li2 molecules with 
zero momentum a) in the BEG limit at 571 Gauss and b) in 
the BGS limit at 892 Gauss of the crossover occurring near 
the broad Feshbach resonance of ^Li at 834 Gauss. The Fermi 
energy of the gas is 380 uK. 



atomic density does, however, not contain the paired 
atoms associated with the dressed molecules that are not 
in the Bose-Einstein condensate. This omission is re- 
paired by the fluctuation corrections which contain both 
the associated changes in the atomic density and twice 
the total density of dressed molecules that are not Bose- 
Einstcin condensed, i.e., twice the dressed molecular de- 
pletion. 

We are now in a position to determine Z . In princi- 
ple it is equal to the residue of the pole of G'ni;ii(0, icj) 
at = 0. To understand the physics of that result bet- 
ter it is useful to consider the spectral function of the 
bare molecules with zero momentum, i.e., pin(0,w) = 
— Im[G'i„;ii(0, a;-|-i0)]/7r. This spectral function is closely 
related to the density of states of the molecules at zero 
momentum and thus gives detailed information on the 
dressed molecular content of the gas. Most importantly 
for our purposed, the Bose-Einstein condensate of dressed 
molecules gives rise to a delta-function in the spectral 
function exactly at zero frequency. The strength of this 
delta function is precisely Z, because this is the probabil- 
ity to take a bare molecule out of the Bose-Einstein con- 
densate of dressed molecules. Besides this bound-state 
contribution, the spectral function contains also a contri- 
bution from the continuum of atomic scattering states. 
In the BEC limit of the crossover the continuum contri- 
bution only occurs at positive frequency and starts at a 
frequency of — ~ 2h/ma? as shown in Fig. la. In 
the BCS limit the continuum contribution occurs both 
at positive and negative frequencies, which start at a fre- 
quency of about ±2g-\/Z?imc/S, respectively, due to the 
gap that exists for the creation of an atomic quasiparticle- 
quasihole pair. This is shown in Fig. lb. The negative 
frequency part of the spectral function is especially im- 
portant, because it determines the depletion of the Bose- 



Einstein condensate of dressed molecules. Physically it 
represents the dressed molecules that are stabilized by 
the Fermi sea !29|. Making use of the above physical 
picture, we finally obtain the desired result. 
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where S,^"^ = (-i)"5"Sy(0, 0)/9w". 

At nonzero momenta the spectral function is similar 
but now contains two delta functions at the frequen- 
cies ±[.i;iji(k), which have the strength 2'|M,n(k)p and 
-Z|u,n(k)|2, respectively, with |um(k)P - |i;m(k)P = 1. 
This shows explicitly how at long wavelengths our theory 
leads to a Bogoliubov-like theory for dressed molecules 
with a wave function renormalization factor Z . More- 
over, in agreement with the Goldstone theorem, the 
quasiparticle dispersion Wm(k) always turns out to be 
linear at long wavelengths. In the following we therefore 
determine also the associated speed of (second) sound 
throughout the BEC-BCS crossover region. 
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FIG. 2: The solid curve shows the probability Z and the 
dashed curve the fraction 2Znmc/n as a function of magnetic 
field. The Fermi energy h^kY/2m of the gas is 380 uK. The 
data points are from the experiment of Partridge et al. [lo|l . 
In the inset the solid line shows the Bose-Einstein condensate 
fraction of dressed molecules 2nmc/?i and the dashed line the 
contribution of the fluctuations to the total atomic density. 
The vertical lines indicate the magnetic fields where fepla] = 1. 

Results and discussion. — In Fig. 2 we show our re- 
sults for Z and compare with the experimental data of 
Partridge et al. general the agreement is satis- 

factory. This is particularly true at relatively low mag- 
netic fields where Z is determined by two-body physics, 
which is exactly incorporated into our theory. At higher 
magnetic fields the theoretical values of Z are somewhat 
higher than the experimental ones. We believe that the 
reasons for this are twofold. First, the experiment is per- 
formed in an optical trap. As a result the experimen- 
tal data involves an appropriate average over the den- 
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sity profile of tlie gas, which lowers the observed value 
of Z. Second, the photodissociation laser used in the 
experiment has a width which is much larger than the 
Fermi energy of the gas. The laser, therefore, has not 
sufficient resolution to probe only the Bose-Einstein con- 
densate of dressed molecules, and probes also the dressed 
molecules which are not Bose-Einstein condensed. This 
second effect should be especially important at high mag- 
netic fields. To disentangle these different effects, how- 
ever, goes beyond the scope of the present Letter and is 
left for further investigation. 

We also show in Fig. 2 the Bose-Einstein conden- 
sate fraction of dressed molecules 2n^c/'n throughout the 
BEC-BCS crossover. In qualitative agreement with the 
poor man's approach of Ref. the latter fraction is 
always substantial below the Feshbach resonance and be- 
comes negligible only sufficiently far above the Feshbach 
resonance when feplaj < 1. This is an important obser- 
vation, because in our theoretical description of the ex- 
periment of Partridge et al. the molecular probe couples 
directly to the dressed molecules, which act as distinct 
entities in the gas since the atom-molecule coupling is 
much larger than the coupling of the probe laser to the 
bare molecules. In this manner it is most easy to under- 
stand the experimental observation that there is initially 
an exponential (one-body) decay of a large part of the to- 
tal atomic density with a rate that is much smaller than 
the bare molecular photodissociation rate. 
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FIG. 3: The speed of sound as a function of magnetic field. 
The Fermi energy of the gas is 380 nK. The dashed curves 
show the BogoUubov result h^/2nan^/m and the speed of 
sound hkp /y/3m = vf/ \/3 of the Anderson-Bogoliubov mode. 



For completeness we give in Fig. 2 also the Bose- 
Einstein condensate fraction of bare molecules 2Znync/n 
and the fluctuation corrections to the total atomic den- 
sity. As expected the fluctuation corrections are very 
important in the crossover region and become small far 
away from the Feshbach resonance, where mean-field the- 
ory applies. From the fluctuations we also extract the 
speed of sound of the gas, which is shown in Fig. 3. Note 
that in the BCS limit the Anderson-Bogoliubov mode 



is recovered. In combination with the presence of the 
sharp peaks in the spectral function in Fig. lb, this shows 
that also the decoupling of the amplitude and phase 
fluctuations of the Bose-Einstein condensate of dressed 
molecules that occurs in the BCS limit is correctly in- 
corporated. We therefore conclude that the RPA-like 
atom-molecule theory developed here gives an excellent 
account of the subtle interplay between two-body and 
many-body physics taking place at the crossover near a 
Feshbach resonance. 
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